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m Abstract 

The singularities of an irrotational magnetic field are lines of electric current. This prop- 
erty derives from the relationship between vector fields and the topology of the underlying 

, , 3-space. When applied to cosmic fields, it allows for a definition of flux tubes and flux ropes 

as cores (in the sense of the physics of defects) of helical singularities. An interesting feature 
is the quantization of the radii, pitches and helicities of force-free flux ropes. In the special 
case when the total electric current vanishes, a force-free rope embedded in a medium devoid 
of magnetic field is nonetheless topologically stable, because it is the core of a singularity 
of the vector potential. This is reminiscent of vortex lines in superconductors. Magnetic 
reconnexion is also briefly discussed in the same framework. As a general rule the concepts 
of the theory of defects developed in condensed matter physics prove useful in raising new 
questions about the physical characteristics of flux ropes and electric currents in the cosmos. 
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I. INTRODUCTION 



One of the still not fully understood features of cosmic magnetic configurations is the 
universal presence of tubes of magnetic flux, i.e. filamentary domains in which the magnetic 
field is concentrated. Parker [1] chapt 10 p. 208 writes: "The concentration of field into 
isolated flux tubes occurs spontaneously, in opposition to the considerable magnetic pressure 
of the concentrated field. The phenomenon challenges our understanding of the basic physics 
of flux tubes". In Reference [1], Parker goes on to discuss a number of physical mechanisms 
to which the stability of the flux tubes might be attributed. However, in a recent review 
[2] , written thirty years after Reference [1] , his view is that the fibril state of solar magnetic 
fields remains unexplained. 

Flux tubes were first observed in the Sun's corona at the occurrence of eclipses, and have 
since been routinely investigated with coronagraphs that occult the solar disk and reveal the 
relatively low intensity photons emitted by the corona. For a review of corona observations, 
see Aschwanden [3] chapt. 1. Flux tubes have also been observed by crossing spacecrafts 
as major structural features of the plasmoids in the Earth's magnetotail [6] chapt. 6. But 
there is some evidence that flux tubes are ubiquitous throughout the Universe on a variety 
of scales, in stars other than the Sun, in galaxies, and in intergalactic space. Thus, one 
expects that the flux tubes in the Sun provide a model for stars where 'starspots' have 
been observed [2] . 'Twisted trunks' visible in some planetary nebulae have been interpreted 
as twisted filaments confined by a magnetic field [4]. Extragalactic filamentary magnetic 
networks in galactic clusters have also been discovered [5]. 

A number of computer simulations have reproduced the major features of these tubes 
but without giving a clear understanding of their origin. It has been suggested that these 
flux concentrations are due to the expulsion of the magnetic field by eddies [7]. This is 
certainly true, but is not the sole reason, we claim. This behavior has been compared to 
the expulsion of the magnetic flux in superconductors of the second species by Mahajan 
& Yoshida [8]; according to their theory, which assumes a completely ionized plasma, flux 
concentration occurs only for a specific range of the parameters that enter the solutions of 
the MHD equations. However, the phenomenon is too universal to be limited to particular 
solutions. 

We propose hereunder a possible explanation that could work in all circumstances. The 
only prerequisite is that the cosmic magnetic field is irrotational throughout most of the re- 
gion of interest, from which assumption we demonstrate that the singularities (the 'defects') 
of such a universal geometry are tubes of electric current (similar to Birkeland currents 
[9]), and consequently, thanks to the Biot— Savart relationship, tubes of high magnetic flux. 
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Thereby, the tube stability is of a topological nature, which is in no way connected to the 
presence of eddies. 

The assumption of an irrotational magnetic field is frequently invoked in the heliophysics 
literature. Aschwanden [3], chapt. 5, stresses that the coronal magnetic field "can be 
quantified in terms of a potential field which characterizes to first order unipolar fields in 
sunspots, dipolar fields in active regions, and is often used to compute the global coronal 
field with a source-surface model". At the same time it is recognized that there are regions 
of concentrated current, in the form of helically twisted loops, or tubes, where the magnetic 
field reaches values in the thousand-gauss range, while outside these regions the field strength 
is typically less than 5% of these values. We shall see that it makes sense to assume that 
the magnetic field is in fact vanishing outside the tubes. 

It is usual to distinguish between flux tubes with or without magnetic helicity H = 
J A ■ BdV. The tubes with no n- vanishing H, the so-called twisted flux tubes or flux ropes, 
whose origin is still disputed (e.g., produced in the convection zone of the sun, or resulting 
from photospheric flows, or . . . ), are attended by coronal mass ejections (CME's) and are 
thought to be related to the formation of prominences [6]. The question of the nature of a 
flux tube is equivalent to the question of the nature of the core of a singularity, i.e., how 
the medium escapes the presence of a physically impossible singular solution which would 
otherwise be unavoidable on topological grounds. The core is a domain with properties 
which render the singularity virtual. Thus, to solve for a singularity is to solve a problem 
where two solutions are in contact at the core boundary. Such a question occurs frequently 
in the theory of defects in condensed matter physics [10]. 

The magnetic-field singularities we discuss here are distinct from those arising from tan- 
gential discontinuities (current sheets), which are a necessary consequence of the equilibrium 
balance of Maxwell stresses [11], as well as the separatrices and null points that shape the 
global magnetic geometry at the surface of the Sun [12]. 

In this paper we shall stress the physical properties of flux tubes and, especially, flux 
ropes in relation to their topological properties. The calculations are developed for cylindrical 
tubes, Section IV, V, VI, and VII; this geometry constitutes a simple example which does not 
require a full understanding of the distinctive analytic and topological features of magnetic- 
field singularities in general. Closed tubes (and their corresponding closed currents) require 
such a deeper understanding. Some general considerations are given in Section III. 

II. TUBES OF ELECTRIC CURRENT ARE TOPOLOGICAL DEFECTS 

Let us assume the existence of a magnetic field B with no associated current (VxB = 0) 
deriving from the following scalar potential: 




(1) 
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The isopotential surfaces are ruled helicoids and the magnetic field, expressed in cylindrical 
coordinates (r, 6, z), is given by 

B = V^ = {0,B o ^,B o }. (2) 
Lixr 

In this note, we wish to investigate the physical nature of the singularity of the field along 
the z-axis. 1 

One encounters a similar helical geometry when investigating screw dislocations in liquid 
crystalline smectic phases. In the ground state, the structure of these phases consists of 
planar equidistant lamellae; they can suffer certain singularities [13] such as line defects 
(dislocations, disclinations), focal conic domains, singular points, etc., whose topology is 
related to the symmetry group of the phase. Screw dislocations are among the simplest of 
these singularities, in which the lamellae take the shape of helicoids. 

Let d be the distance between lamellae in the ground state of a smectic phase. It follows 
that the pitch p of the helicoids, which are given by ipi am — z T (nd/2n)6 = 0, n G Z, is a 
multiple p = ±nd of this distance. The pitch can be expressed as an integral along a loop 
C (the Burgers' circuit) surrounding the z-axis (the locus of the dislocation), namely 

j> u ■ ds = Aip lam = P, 

where u = V^ am = {0, ^fnd/2irr), 1}; pz is then the Burgers' vector. 

Analogously, consider the integral § C ~B ■ ds along a circle C of radius r. One gets: 

B • ds = B b, (3) 

c 

which is independent of the circuit C surrounding the z-axis. This quantity measures the 
current intensity X through C, namely (An/c)! = /JVxB • da; VxB vanishes almost 
everywhere, so we use the language of distributions, 

4-7T 

VxB = —XS(x)5(y) z, (4) 
and, by employing Stokes' theorem: 

47T 

-I = B b. (5) 

c 

The current intensity appears as the equivalent of the Burgers' vector. It is concentrated 
along the z-axis. But by spreading the current in a cylindrical region of finite radius r , 
we get rid of the singularity, and the problem now is to find an irrotational magnetic field 
outside the cylinder and a rotational magnetic field inside (the core region) subject to the 
boundary conditions [B ■ n] = and [A] = (see [14] chapt. 5) and the global condition 

pro 4ft 

2tt / rdr z • VxB = — X = B b. (6) 
Jo c 



1 Certain developments hereunder are inspired by the theory of defects in Condensed Matter Physics; the 
reader will find a short reminder of the main concepts of this theory in the appendix A. 
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Other conditions may apply in the resolution of the fundamental equations of magnetohy- 
drodynamics, but we shall restrict here to steady-state solutions with no material current, 
i.e. to a purely magnetostatic problem. We also assume that the medium is a perfect con- 
ductor. Finally, we impose cylindrical symmetry, so that B depends only on r. 

Remark 1: The current along the z-axis is of course the source of a magnetic field B 
through the Biot— Savart relation. 2 This is not the way we have introduced it; we have 
instead introduced a set of surfaces orthogonal to an irrotational magnetic field. At this 
point, we wish to stress that such a set of surfaces cannot be defined unless the magnetic 
field has nonzero curl; a sufficient condition for the existence of a set of surfaces orthogonal 
to a given magnetic field B field is precisely that B be irrotational; in effect, VxB = => 
v ■ Vxf = 0, where B = |B|za 

Remark 2: The preceding results are topological in essence. The constancy of nd in the 
lamellar case and of X in the electric-current case are physical properties that have to be 
maintained whatever the deformation of the helicoidal surfaces with respect to the regular 
helicoids may be, provided the helicoidal topology is preserved. In the electric-current case, 
X being given, one can smoothly bring the actual potential surfaces, in the vicinity of a 
straight line current, to the shape of regular helicoids (which ensures h = est and Bq = est), 
such that (5) is satisfied. 

Remark 3: The analogy between the topological defects in smectics and the topological 
defects of an irrotational magnetic field is restricted to screw dislocations, but merits discus- 
sion because it yields the main characteristics of line defects in magnetic fields, in particular 
their helical features. Furthermore, magnetic configurations in the cosmos can exhibit cur- 
rent sheets [11], which have no analogs in smectics (such singularities cannot be foreseen 
from the symmetry group of the lamellar phase; the normal v to the lamellae is not a vector 
(nor a pseudo-vector, like B), but a director, i.e., v = —v, see de Gennes [15]). Magnetic 
defect configurations in ferromagnetic materials, with domain walls and Bloch and Neel 
lines, are somehow topologically analogous, but with nothing comparable to the presence of 
electric currents [16]. On the other hand, vortex lines (vortex sheets) in irrotational fluids 
are topologically analogous to current intensity lines (current sheets). 

III. THE NATURE OF THE ELECTRIC CURRENT SINGULARITY 

The analysis which follows is an application of the relationship between the topology of 
three-dimensional domains and the calculus of vector fields defined on such domains (see 
ref. [17] for a review). Our starting point is (a special case of) the Hodge decomposition 

2 More precisely, the azimuthal component of the magnetic field is generated by the z-component of the 
current, and the azimuthal component of the current contributes to B z inside the flux tube. But the 
(constant) z-component of B outside the flux tube is not generated by the current, but is generated by 
external currents that are far away. We shall later consider a situation with no magnetic field outside the 
tube of current, in which case the Biot— Savart relation completely determines the field inside the tube. 



theorem. Let u be a (sufficiently smooth) vector field on a domain C 1R 3 which satisfies 
tangent boundary conditions, i.e., u n = 0. (Here, u„ denotes the normal component of u on 
the boundary dfl, and u t below denotes the tangential component.) Then u can be uniquely 
expressed in the form 

u = Vxa + Vp + h, (7) 

where a has zero divergence and satisfies normal boundary conditions while its curl satisfies 
tangent boundary conditions, i.e., 

V-a = 0, a t = 0, (Vxa)„ = 0, (8) 

Vp satisfies tangent boundary conditions, i.e., 

(Vp) n = 0, (9) 

and h is a harmonic vector field satisfying tangent boundary conditions, i.e., 

V-h = 0, Vxh = 0, h n = 0. (10) 

The conditions on h imply that Ah = 0, where A is the Laplacian. Moreover, a and h are 
uniquely determined by u, while p is defined up to an additive constant. The components 
of the Hodge decomposition (7) are mutually orthogonal, i.e., 

/ (Vxa) • Vp = / (Vxa) • h = / (Vp) • h = 0. (11) 
Jn Jn Jn 

We note that the divergence and curl of u are given by 

V • u = Ap, V x u = -Aa. (12) 

Thus, the Hodge decomposition implies that if B is an irrotational magnetic field on f2 
satisfying tangent boundary conditions, then B is harmonic. 

It turns out that harmonic vector fields satisfying tangent boundary conditions are closely 
related to the topology of the domain on which they are defined. More precisely, the space 
of such fields is finite dimensional, with dimension given by the genus g of the boundary dfl 
(or, if dQ consists of a number of components, the sum of the genera of its components). 
Indeed, harmonic vector fields can be characterised by a set of g invariants, either fluxes or 
line integrals, as follows: Construct a minimal set of surfaces Ej in Q such that Q — UEj is 
simply connected - in the language of defects, these are 'cut surfaces'; mathematically, they 
are Seifert surfaces. The number of cut surfaces is given by g. The intersections of the E 3 -'s 
with dQ yield g closed curves, d'Ej, on the boundary dQ. One can also define g additional 
closed curves, C k , on dil such that C k intersects <9E fc once and does not intersect any of the 
other dEj's. Together, the dE^'s and C^'s constitute a basis for the first homology group of 
the boundary dfl, which therefore has dimension 2g. 
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Then, and this is the basic point which allows one to make contact with the theory of 
defects, a harmonic magnetic field B is uniquely determined by its line integrals around the 
Cfc's (which may be interpreted in terms of currents outside Q), or, equivalently, by its fluxes 
through the cut surfaces £j. As explained below, there are physical reasons to emphasise 
the first determination, namely 



In simple examples this construction is easily visualised. For example, Mahajan & al. 
in ref. [18], for the purpose of modelling a tokamak, take Q to be a solid torus T; S, a 
cross-section of the torus; <9£, a closed curve going round the torus the "short way"; and C, 
a closed curve going around the torus the "long way". In this case, g = 1. 

The main example we consider, also with g = 1, is directly related to the notion of defect 
developed in this paper, and is shown in Figure 1. The domain Q is obtained by removing 
a solid torus T from the interior of a (deformed) three-dimensional ball. The boundary 
dfl consists of a simply-connected outer component, dfl out , which is just the boundary of 
the three-dimensional ball, and an inner component, dVL in = ST. In order to render Q 
connected, we introduce a cut surface £ whose boundary <9£ is a closed curve which goes 
around T the "long way". For ease of notation we will denote <9£ by L. We also introduce 
a closed curve C which goes around T the "short way", which intersects L once. 



Figure 1. A domain Q with g = 1. The colored torus is supposed to be empty (this restriction can 
be lifted; then the inside of the torus is the core of the singularity). dT, is along the 'long way', 
and C along the 'short way'. The cut surface X is shadowed. 

The quantity (3 given by (13) takes the same value if C is replaced by a homologous loop 
C since, by Stokes theorem, the difference § c B • ds — § c , B • ds is given by the integral of 
V x B over a surface in Q, and V x B vanishes by assumption. This invariance gives C 




(13) 
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the status of a Burgers' circuit, and dT the status of a line defect, inasmuch as the radius 
of T can be reduced to zero in such a way that T collapses onto the closed curve L. In this 
limit, B becomes singular on L, and L has the characteristics of a vortex line, with E as cut 
surface. 

We may regard (3 as proportional to a current X threading the interior of the torus T. 
Indeed, in the limit that T collapses to L, we may write that 

VxB= ^XL5f(r), (14) 

where L is unit- vector along L and <5|(r) denotes the normalised two-dimensional delta 
function on L. 

As B is irrotational, it may be expressed (locally) as the gradient of a (multi-valued) 
scalar field ip. Introducing local cylindrical coordinates (r, 9, z) around L, we may write 

ip = —6 + 72 + 0(r 2 ) + higher order terms. (15) 

There is no linear term in r, since B„ = on dT. (1) is a particular instance of (15), where 
Q is the solid cylinder with its axis removed (and the higher-order terms in (15) are absent). 

Finally, we note that provided the current X remains fixed, (3 is invariant under deforma- 
tions of L. This is in contrast to the alternative invariant, namely the flux of B through the 
cut-surface S, which depends not only on X but also on the geometry of Q. As current-tubes 
in the cosmos are expected to be mobile while keeping their current constant, we choose (3 
as the characteristic invariant of a harmonic magnetic field B. 

In the following sections, where we restrict to straight tubes of flux, the full characteristics 
of this general theory of singularities in an irrotational magnetic field will not be necessary, 
and it will be sufficient to think of the tubes as the cores of helical defects with regular 
helical level sets ip. 



IV. THE ISOLATED FLUX TUBE 

As a first, instructive example, we briefly consider what is possibly the simplest model 
of a cylindrical defect, in which the current density (4ir/c) zj = z- VxBis uniform in the 
core, i.e. (Att/c)z ■ j = B^)ji\r\^ where r is the core radius. As we shall see, this model is 
too simple in several respects. First, it is not a static solution of the MHD equations (unlike 
the examples in later sections, this configuration is not force- free); second, it leads to an 
unphysical dependence of the core radius on field strength. 

Since the current density vanishes outside the core, there may be a surface current 
(47r/c)j >5 . = [n x B], but we first assume no tangential discontinuity of B at the core bound- 
ary. Thus B is continuous. (2) gives the field outside the core, and the conditions just 
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indicated yield the field inside: 

The total magnetic energy per unit length for a cylinder of radius R > r is given by 
With B as above, we obtain 

„ B 2 ( b 2 b 2 R 2 \ 

£ = - 1 — + — In — + ttR 2 ) 18 

In the absence of additional terms in the free energy, the energy decreases without limit 
when tq — > R. That is, the core occupies all available space. This example is therefore of 
very limited interest. Furthermore, even for r < R there is no flux concentration in the 
core; indeed, the flux per unit area is just B both inside and outside the flux tube. 

Let us consider a less restrictive model, which includes Parker's model and the one above 
as particular cases. We allow for a nonvanishing field inside the flux tube of the form 

B = Bo f 0,-^,1^ (19) 



27rr, 



The field outside is still of the form given by (2), but of strength fii different from fi , so 
that B may be discontinuous across the core boundary, with corresponding surface current 
density given by 

J S =±(B, -*)(«- A..) (20) 

Flux concentration corresponds to \B \ > \Bi\. 
In this case, the energy is given by 

„ B 2 (b\ 2 B 2 (b\\ R B 2 -B 2 B 2 2 , , 

£ = — [ — ] + — — n — + ^r 2 + -^-R 2 . 21 

16 \2irJ 8 \2irJ r 8 8 V 

As a function of the flux-tube radius, the energy displays a minimum (d 2 £/dr 2 > 0) for 



^0 = 



\B 2 -B 2 2n> 
The minimum energy can then be written as 



Bl 161 if |Bb|>|Bi|. (22) 




Thus, the assumption of flux concentration is compatible with an energetically stable con- 
figuration. If I S 1 3> I one gets 

-!H(t) 2 )!- 

The helicity vanishes identically inside the core, and vanishes also outside in the limit fii = 0. 
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V. THE TWISTED FLUX ROPE 



Static configurations of currents and fields may be obtained by requiring the Lorentz 
force to vanish, i.e, j x B = 0. Hence we consider cores with force-free configurations. This 
is a particular model of a flux rope in the sense of Parker [1] chapt. 9, although Parker 
assumes that there is no magnetic field outside the tube; we consider this case, which yields 
the most physically sensible results, in Section VII. The force-free field situation is plausible 
for the magnetic ropes of the corona, where the (3 number is very small (/3 is the ratio of the 
plasma gas pressure, n e kBT, to the magnetic pressure, _B 2 /87r), so that the pressure gradient 
is negligible compared to the Lorentz force, which therefore must vanish in a magnetostatic 
situation. The condition for a force-free field is given by the Beltrami equation, 

IVxB = B. (25) 

For the case where i is taken to be constant, the solution in a cylindrical domain is well 
known since Lundquist (see Alfven & al. [19]); it can be written in terms of Bessel functions 
of the first kind: 

Be = AJ X (0 , B z = AJ . (26) 
Requiring B to be continuous on the boundary r = r , we get that 

AJ 1 (r,) = ^, AJ (r)) = B , (27) 

where rj = r Q /£, which implies that 



J (r)) 2-nr ' 



(28) 



We assume that the two essential invariants of the flux rope, namely the current intensity 
X = (c/4ir)B b and the magnetic field B , are given data; i appears as a function of r 
through the boundary condition (28), and is an unknown of the problem, which is to be 
determined, as is ro, by minimizing the energy. In what follows, without any essential loss 
of generality, we may restrict attention to t > 0, while b may be of either sign. 

Notice that, when navigating around the cylinder one turn, i.e., a distance 27rro, a repre- 
sentative point at z = zq is lifted to z + b along a geodesic of the cylinder, see Figure 2. Thus 
tan0 = Ji(r]) / J (r]) , where tan0 is the slope of this geodesic, which develops into a helix 
with pitch —b on the cylinder r = r . Therefore the lines of force of the magnetic field, which 
are orthogonal to this helix, are helices of opposite sign. If e.g., the helicoids are left-handed 
(b > 0), the magnetic field follows right-handed helices of pitch p(r ) = (27rr ) 2 /6 at r = r . 
And for any r < r Q , the pitch of the helical lines of force of the magnetic field is given by 
p{r) = 2nrJ Q (r/£)/Ji(r/£). Thus, depending on the value of r/£, the pitch may alternate 
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Figure 2. The helical intersection of an helicoid ijj = est with the boundary of the tube, b > 0. Development 
on the plane. 

in sign inside the cylinder. Indeed, if we let ao,j, j = 1, 2, 3, . . ., denote the zeros of Jo, and 
ax tJ similarly denote the zeros of J\, then p(r) is positive on the intervals &o j3 < r/i < aij 
and negative on the intervals a 0j+ i < r/i < ai i3+ i (note that the zeros of J and J\ are 
interlaced, i.e. a 0j < ai i3 < a 0j +i, and that the magnetic field is horizontal at r/i — a 0j and 
vertical at r/i — aij). See Figures 3, 4, and 5 hereunder. 

The number of sign changes of the pitch is even or odd according to the relative sign of 
the interior and exterior fluxes. Indeed, outside the tube, the magnetic flux per unit area 
is given by Bq. The flux inside the core, $j = JJ^ B(r)rdrd6, follows from the force-free 
condition and Stokes' theorem, 

= MB . (29) 

The interior and exterior fluxes have the same sign if b > and opposite signs if b < 0. At 
the centre of the core, the pitch of the magnetic lines of force is necessarily positive (since 
£ > 0). For b > 0, the pitch of the magnetic lines of force is positive at the boundary; 
therefore, if the the pitch changes sign in the core, it does so an even number of times. 
For b < 0, the pitch of the magnetic lines of force is negative at the boundary; it therefore 
necessarily changes sign inside the core (an odd number of times). 
We remark that i can be given a simple interpretation. Because 

11 

- = jgpB • VxB = v ■ Vxi/, where B = f|B|, 

w = —2tt£ is the pitch of the vector field u, which thus rotates helically with a constant 
velocity along the radii of the cylindrical geometry. One can say that B suffers a double- 
twist, one along z, the other along r, which is analogous to the double-twist observed in 
liquid crystalline blue phases and cholesteric phases, where in general w depends on r. Of 
course w can also depend on r in the MHD case, and in a more general geometry probably 
does, as a rule, without breaking the force-free condition. 
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VI. MAGNETIC ENERGY OF THE TWISTED TUBE 



The magnetostatic energy (17) is the sum of the exterior contribution 



Bn ( 1 , „q o, ( b \ , i? 



and core contribution 



^=^\\{R 2 -rl) + \- \ ln^| (30) 



where we have used the identities [20] 

J\dpJl(p) = \f (4(n) + Jl(v)] 

jf%4> JfM = \<f [m + J?M - 2^M) . (32) 
We obtain 

„ 5 2 / 2 2 6 2 , i? 6 2 &A 

£ = — LR 2 + r 2 + — In — + — , 33 

We assume that b and So are independent invariants, and seek to minimize the total 
energy £ with respect to tq. Recall from (28) that the Beltrami length i depends on ro; for 
convenience, we write the relation in the form 

F = r °f {ro/e) -2n = °- (34) 

By the chain rule, one has 

d£/dr = -(dF/dr )/(dF/d£). (35) 

Noting that 

fw=-|(!)V), Iw=(|)(-/^Kf)'^»' < 36 > 

and that 

one obtains, after replacing J\J Jq by (cf (34)) 



C = ^~ (38) 
2vTr n 



the following relation: 

dt c + r 1 



dr l-( ro /£)(C + C -i) 
12 



(39) 



The extrema of £ are given by dS/dro = 0, i.e., 

dr 4 °\ ^H-( ro /£)(( + (-i) J ~ 
After some easy algebraic manipulation, this equation can be written: 

b 2C 2 



0. 



(40) 



2tt£ 1 - C 4 



(41) 



It follows that if b > 0, then < ( < 1, while if b < 0, then ( < — 1. For ( = ±1 one finds 
that £ = 0: there is no twist. 




it) n = r / >■ 



Figure 3. Nonvanishing magnetic field outside the rope and continuous at the boundary, with £ > 0. Blue 
curves: boundary conditions at r — ro; other curves: plot of the minimization condition d£/dr$ = 0, 
C > for b > 0, £ < for b < 0. See Tables I and II for numerical values. 



The relation between ro and I can be better understood by consideration of Figure 3. 
The blue curves are plots of the boundary condition ( = Ji(j]) / Jo(rj) (cf (34)), and are 
asymptotic to the vertical lines i] = ao, r The other curves are plots of the minimisation 
condition (41), expressed in the equivalent form 



n 



2( 



i-C 4 



(42) 



The upper red curve corresponds to b > 0. It is tangent to the graph of Ji(v) / Joiv) a ^ 
i] = 0, and has asymptote ( = 1. Thus, for rj large, one gets solutions r ~ b/2n, whatever 
the value of £ may be (see Table I). The lower red curve has asymptote ( = — 1. Thus, for r] 
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Figure 4. Nonvanishing magnetic field outside the rope and continuous at the boundary, with t < 0. Notice 
the change of the sign of b with respect to Figure 3. 

large, one has solutions , r ~ —b/2n (see Table II). The case £ < is illustrated Figure 4; 
the symmetries with respect to Figure 3 are obvious. 

Thus, the critical values of l 3 and tq 3 are quantized for Bq and b given. For b positive, 
T} 3 = r^j/ij lies in the interval (aij,aoj+i); for b negative, in the interval (aoj,ai i: ,). The 
second derivative d 2 S/dr 2 is positive all along the curves d£/dr = 0. 3 Thus the solutions 
are minimizers of the energy. The first interval is special. The origin is indeed a minimum, 
but with r = oo. 

As noted above, for large |t7j|, we have that ro,j ~ \b\/27i. Thus, for large j, the energy 
(33) is given approximately by 

Bib 2 ( 2 , 27riA B 2 

£ 7 ~— 1 + 2 In—— +— -R- 43 

3 8 4vr 2 V V 3 \ b \ ) 8 

It is straightforward to obtain the following asymptotic expressions for r) 3 , 

0-1/2)*- a^jri), j>i,6>o, 

rj 3 w jk - — , j > 1, b < 0, (44) 



3 One gets, starting from eq. (40), that d£/dr = (Bg/4)r (l - C 2 + C(C + C _1 )/-D) = 0, where D = 
1 — (r /£)(C + C" 1 )- Using (41), one gets that, at a critical point r of £ , 

d^ = _£| 2 2 ( C + r i) 
9r 2 4 £>r ^ ^ j 

Using (41), it is easy to show that bD < 0, regardless of the sign of b. Therefore d 2 £/dro 2 > at all 
critical points. 
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from the asymptotic expressions for Bessel functions of large argument [20] , 



Jv 



,{x) = •sjlj (iTx) COS — -VII — -TT^j + 0(l/\x\), \x\ ^> 1. 

The limiting value of the energy £ 3 as j — > oo is given by 



(45) 



8 4tt 2 



1 + 2 In 



2ttR\ B, 



D 2 



\b\ 



(46) 



Tables I and II present some numerical values attached to the flux tubes, both for I > 0. 
The last column relates to the magnetostatic energy, except for the logarithmic term, which 
varies more slowly than the others, and the R 2 term, where R may be understood as a mean 
distance between neighbouring ropes. The fact that this last term is so large might explain 
the preferred occurrence of the model of the next section, in which this term is absent the 
free energy. 



J 




c- - b 




( 2 f ? (r 2 o ~ b i) 


'h — i 




b 


2 


4.6124 


0.8859 


0.2447 


0.7848 


3 


7.7995 


0.9338 


0.1373 


0.8721 


4 


10.9581 


0.9533 


0.0957 


0.9088 


5 


14.1087 


0.9639 


0.0735 


0.9292 


10 


29.8321 


0.9831 


0.0341 


0.9665 


20 


61.2548 


0.9918 


0.0165 


0.9837 


30 


92.6729 


0.9946 


0.0108 


0.9892 


50 


155.5064 


0.9968 


0.0065 


0.9936 


100 


312.5873 


0.9984 


0.0032 


0.9968 



Table I. £ > 0, b > 0; the term in In R/b is omitted — R is an unknown parameter, and the In term 
is varying more slowly than the other terms 



The magnetic flux is Bq per unit area outside the core, and B^hijur^ = Bq(1 — £ 4 ) inside 
the core. There is therefore no concentration of the magnetic flux. The largest positive 
value of 1 — C 4 , namely 0.384 obtains for £ > 0, b > 0, j = 2 (Table I). For this solution, 
the magnetic lines of force change handedness once, which of course decreases the flux total 
through the rope. The largest-in-magnitude negative value of 1 — £ 4 , namely —0.754, obtains 
for £ > 0, b < 0, j = 2 (Table II). This negative density of flux inside the rope, which is of 
opposite sign to the outside flux, exceeds in modulus that of the preceding case b > 0. In 
any case, the flux concentration is very small, and even absent for j > 2. 
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3 


„ _ !Xj 
VJ — £ 


c- - b 


2nl, 
b 


( 2 f ? (rg - b i) 


1 


3.0526 


-1.1508 


-0.285 


1.3247 


2 


6.2320 


-1.0770 


-0.1490 


1.1600 


3 


9.3889 


-1.052 


-0.1013 


1.1064 


4 


12.5388 


-1.0391 


-0.0768 


1.0797 


5 


15.6855 


-1.0314 


-0.0618 


1.0637 


10 


31.4044 


-1.0158 


-0.0314 


1.0318 


20 


62.8260 


-1.0079 


-0.0158 


1.0159 


30 


94.2438 


-1.0053 


-0.0106 


1.0106 


50 


157.0772 


-1.0032 


-0.0063 


1.0064 


100 


314.1581 


-1.0016 


-0.0032 


1.0032 



Table II. I > 0, b < 0; the term in lni?/|6| is omitted. 



VII. THE TWISTED ROPE WITH NO FIELD OUTSIDE 



We consider an isolated cylindrical magnetic rope of radius r with force-free internal 
field and vanishing external field. The internal field is of the form 

B e (r) = AJ^r/i), B z (r) = AJ (r/£), (47) 

where £ is the (constant) Beltrami parameter. It is convenient to express the parameters 
A and £ in terms of physical parameters, namely B , the flux per unit area just inside the 
boundary, and X, the current inside the rope. The relations are then the same as those 
obtained in the preceding sections, namely 

A = B J (r /£), ^( ro /*) = -*-, (48) 
J 2nr 

where 

The discontinuity in the field across the boundary requires a surface current density, which 
is given by 

* = h = -^- (50) 

This model is commonplace in standard treatments of magnetic ropes. The new feature 
here, as in the preceding section, is an additional constraint which derives from an underlying 
helicoidal topology. In addition to the interior current X, we regard the pitch of the field- 
line helices as an invariant. Note that, as the external field vanishes, the helicoidal topology 
is virtual, and should be understood as arising from the limit of field configurations in 
which the ratio of the external to internal field tends to zero at constant pitch (so that the 
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helical character is preserved). Alternatively, the helical structure may be derived from the 
vector potential, as follows. Inside the core, the vector potential is £B (c.f. (25)), which 
is divergenceless. Therefore, by continuity, the vector potential outside the core is given 
by A = £(0, B()b/(27rr), Bq). 4 As the external vector potential is irrotational, it may be 
expressed, locally at least, as a gradient V0, with = B £(z + (b/2n)9). The level sets of 
<f> are then helicoids of pitch b. In this way, the twisted rope may be viewed as the core (in 
the sense of defect theory) of a helicoidal topology with vanishing external field. (Note that 
under a change of gauge A — > A + V/, the scalar field <fi becomes <fi + /. The helicoidal 
level sets are deformed, but as / is necessarily regular (A is required to be continuous), their 
topological structure, and in particular the core, persists.) Taking B and b to be invariants 
of the magnetic rope is then equivalent to taking the current and pitch as invariants. 
The energy in this case is calculated similarly to (33), and is given by 

The extrema are characterised by the condition dS/dr = 0. With d£/dr given by (39), 
this may be written as 

< 4 +M) c2 -^°' (52) 

V = + - *V 2 , where rj = r /£. (53) 



or equivalently, 



2 s 1 + C 2 ' 
From (48), we have the consistency condition 

C=yfa)- (54) 

Critical points are then given by simultaneous solutions of (53) and (48). These are shown in 
Figure 5, and numerical values are presented in Table III. In each interval &i, 3 -i < r /£ < a, Qj , 
there is a single critical point r 0j with corresponding Beltrami parameter £ r 

In fact, all critical points of £ are local minima. A straightforward calculation yields, 
after some algebra, the following expression for the second derivative of the energy evaluated 
at the critical points: 

d 2 £ _ B 2 C 2 - * 



dr 2 2 C 2 + l' 
It is readily verified from Table III that ( 2 > 3 for all critical points. 



(55) 



4 The vector potential being continuous amounts to saying that there are no (^-function singularities in the 
B-field. Consider e.g. the discontinuous vector potential 

A(r,9,z) = 0, r < r , A(r, 9, z) = CV9, r > r , 

then B = Vx A = C /{2irr )8{r — r )z. That is, one gets a (5-function B-ficld on the cylinder r = r . 
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For large j, rjj is large and Q ~ 2r]j. From the asymptotic expression for aoj (which 
follows from (45) and the fact that Jg = — Ji), one readily obtains 

b „ b 



'^- 47r2 y ^~ 47r3 yr ( 56 ) 

Thus, for large j, the radius of the flux rope goes to zero even as the number of oscillations 
in the handedness of the field increases. The energy (51) in this limit is given by 



t = J 1 /J [=b/Urr | 



J 




10 

5 

a l,D 








-6 




r. 

-10 


2 




6 

1 



Figure 5. No magnetic field outside the rope; the two first solutions for I > and t < 0, corresponding 
respectively to 6 > (right-handed magnetic field lines of force) and b < (left-handed lines of force). See 
Table III for numerical values. 

E ~ (l - ~ ^ 1 

J ~ 16tt 2 V 4ag ;J ) / ~ 16vr 2 ^ 4vr 2 j 2 

The energy scales expected, and increases but slightly when j increases, and thus 

does not depend strongly on the radius of the rope. The energy of the flux rope with 
vanishing external field is clearly smaller than the case considered in the preceding section, 
where the surface current was taken to vanish and B approached the uniform field B z, 
even if the energy of the uniform field Bq z is neglected. 

The magnetic flux density in the core is B b£/irrl = 4_B (1 + C 2 )/(3 + C 2 )- Thus, for 
a given amplitude B , the flux density tends very quickly to 4B Q for £ — y oo; in fact it is 
already ~ 3.53-Bo for the first quantized state j = 1. In this state the magnetic lines of force 
do not change handedness in the core. But the flux concentration increases even when j 
increases, i.e. for states whose handedness changes sign 2(j — 1) times. 

In the preceding calculation, the twisted rope has been treated as the limit of a rope with 
B — y with fixed interior current. In a forthcoming publication this assumption will be 
discussed and compared to an alternative, namely a flux rope with fixed internal magnetic 
flux. 
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J 




c- - b 

27rro l7 


2n£, 
b 


1 _)_ ( 27rr O,j 2it£j 


1 


2.1288 


3.7610 


0.1249 


0.9458 


2 


5.4258 


10.6657 


0.0173 


0.9915 


3 


8.5950 


17.0733 


0.0068 


0.9966 


4 


11.7488 


23.4123 


0.0036 


0.9981 


5 


14.8973 


29.7273 


0.0023 


0.9989 


10 


30.6183 


61.2039 


0.00053 


0.99973 


20 


62.0404 


124.0647 


0.00013 


0.99994 


30 


93.4584 


186.9060 


5.725xl0~ 5 


0.99997 


50 


156.2918 


312.5773 


2.047xl0~ 5 


0.999989 


100 


313.3727 


626.7421 


4.0916xl0~ 6 


0.999997 



Table III. The last column is the line energy, to a factor j^B^b 2 . 

VIII. DISCUSSION 

The foregoing gives a very idealised picture of a flux tube, even if it correctly captures the 
singular nature of a tube of current, i.e. that the current is a singularity of the magnetic field, 
or, in the case of a twisted rope with no field outside, an irrotational vector potential. And 
indeed, the true physical invariant is the current intensity Bob, which stems from the no- 
divergence property of the current density in the static case. However, the results obtained 
are of some importance in the analysis of the physics of flux ropes, in particular those where 
the magnetic field is confined (no magnetic field outside the core, Section VII), to which 
we restrict attention in what follows. 5 A number of solar observations of coronal loops 
seem to indicate that segments of constant cross-sectional radius, implying constant b, can 
exist along long distances; see figures in [22]. Also, according to [23], flux ropes relax after 
helicity injection to a state of axially invariant cross-section. Our model is consistent with 
these observations. 

The equations we obtain are based on a minimization process at constant current, eq. (5) 

X = (c/A7r)B b. 

But a flux rope also carries two other invariants, namely the magnetic flux eq. (29) 

$ = B b£, 

and the helicity H per unit length. This latter quantity can be written 

H = J A ■ BrdrdO = 8ir££, (58) 

5 Remember that a magnetic vortex line in a type II superconductor is also the singularity of an irrotational 
vector potential; the magnetic field is expelled from the specimen, except in the core region that carries 
a quantized magnetic flux [21]. 
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as we may take A = £B for a force-free field. These three invariants are not linearly 
related, whereas they add linearly when several ropes gather to form a single rope. Therefore 
any phenomenon implying a splitting or a reconnection of ropes is necessarily attended by 
dynamical processes, which result in the modification of some or all the invariants. We 
briefly discuss a few examples, where two of the three invariants are conserved during such 
dynamical processes as reconnection or splitting. 

1— X, TL conserved. We assume that a j-rope, with j large, current intensity (c/ArfB^bj, 
is split into j 1-ropes belonging to the j — 1 state, each with current intensity (c/47r)5( 1 )6 1 . 
The conservation conditions can be written: 



Notice that, according to Table III, the energy per unit length is not very different from one 
state to the other, for j large enough; hence the coefficient a is nearly unity. Comparison of 
the two conservation laws yields 



so that the radius of the subdividing ropes is not much different from the radius of the 
original rope. Finally, assuming B^ = -B(i), the Burgers' vector of the subdividing ropes is 
j times smaller than the Burgers' vector of the original rope. 

The energies are considerably reduced in this process; the energy £ after splitting is: 



which favors the process, but at the expense of a large increase of the flux; one finds 



jB{x)bi = B {j) b 3 , £ 3 {B b) b 3 ) 2 = aMB^hf, a « 0.9458 



(59) 



One finds, using (44) for r]j and Table III for rji, that 



ri « 0.7164 r 3 



£ = ] £i = (a/j) £ : 



& = ]&i = (j/a) $ 



2— $, X conserved. The conservation conditions can be written: 



(60) 



These equations yield 




and a decrease of energy 



£ = 3 £i = (a/j) £j- 



This process is thus feasible, at the expense of a change of helicity 
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'H = jH 1 = (a/j) H 



Notice that, using (44) for r\ 3 and Table III for 771, one finds 

r 1 ps (0.6965/j) r r 

3— %, $ conserved. This is the only case out of the three considered here where the 
splitting is not favored; on the contrary, one expects that the elementary ropes reconnect 
to form a larger one, with a large dissipation of energy. This possibility might be related 
to a question raised by Parker [2], namely "what causes the subsurface convection of the 
Sun to sweep thousands of magnetic fibrils together and then to compress them into two or 
more kilogauss?" Note further that the conservation of helicity is an assumption often made 
in high electrical conductivity media [24], while the conservation of flux is also physically 
reasonable. 

Start from a j-rope split into j 1-ropes. The conservation conditions can be written: 

jAB(i)6i = tjBtfbj, e 3 (B {3) b 3 f = aMB^h) 2 , a ps 0.9458 (61) 
The comparison of the two conservation laws yields 

h = {a/j)l 3 . 

This is the inverse of the relation in the previous case, and one finds indeed that the energy 
is considerably increased in this process 

S = jS l = (j/a)£j, 

which is thus not favored. On the other hand, the gathering of 1-ropes into large j-ropes is 
energetically favoured, in a process in which the current intensity varies, which means the 
appearance of transient displacement currents, i.e. transient electric fields. Using (44) for 
T]j and Table III for rji, one gets: 

r 3 ps 1.5603 j 2 ri. 

These are of course ideal processes, which demand to be explored theoretically for their 
own sake. A real process would certainly be a mixture of these, not to mention the well- 
known kink processes that produce writhing in the flux ropes. 

A number of numeric simulations of single-fluid MHD have been performed, all using 
models of twisted flux ropes with force-free fields and conservation of helicity, e.g. [25, 26], in 
order to better understand the phenomena of reconnection (not splitting) which yield flares, 
prominences, etc, in the Sun's corona in relation with the huge energy which can be released 
in these processes. Various modes of reconnection have been recognized, depending on the 
relative signs of the twist, the relative directions of the mean field b, etc. In this discussion we 
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have limited ourselves to parallel tubes of the same twist. It would be interesting to extend 
our calculations to other cases. It would also be interesting to compare the quantitative 
features of the flux ropes obtained from the foregoing analysis - for example, the structure 
of the interior and surface currents, and their relation to the radius of the flux tube - with 
the structures that emerge from numerical simulations in the high-field regime. 
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Appendix A: Some fundamentals of defect theory 

The key ingredient of the theory of line defects in an ordered medium Q (a crystal) is 
the Volterra process [10]; it describes the creation of a line singularity characterized by an 
isometry of the medium, namely a translational symmetry b and/or a rotational symmetry 
(u) = uju), oj being the angle of rotation, v its axis. 

Consider in a perfectly ordered Q a line L, and cut the crystal along an element of surface 
S bound by L (L = <9£ is a loop or an infinite line); £1 and £2 are the two lips of E, the cut 
surface. Displace the two lips Si and £2 by a relative rigid displacement b and a rotation 
uj . viz. by a total displacement (we assume u small): 

d(r)=b + wAr, (Al) 

b is the Burgers' vector of the line defect, uj its rotation vector [27]. Introduce now a perfect 
wedge of matter in the void left by the displacement, in order to fill it with exact contacts 
with the lips (which is possible since d(r) is an isometry) — or remove matter in the regions 
of double covering. Then reestablish the bonds at the lips Si and £2 and let the medium 
relax elastically. The line defect thus obtained is called a dislocation if oj = 0, a disclination 
if b = 0. Figure 6 shows dislocations in a Smectic phase, in which there is only one finite 
repeat distance do perpendicular to the layers (here |b| = 2d). A screw dislocation is when b 
is along the line L (here along the axis of the helically transformed layers); an edge dislocation 
when b is orthogonal to L — here b is in the plane of the drawing, L perpendicular to it, 
Ibl = d. 




a) 



b) 



Figure 6. Elementary dislocations (displacement d(r) reduced to a pure translation) in a Smectic phase: 
a) screw dislocation, b is in the direction of L; b) edge dislocation, b is perpendicular to L. In both cases 
L is an infinite straight line. 



Such a process has of course no meaning along L itself, and one assumes henceforth that 
a toric region has been removed along this line. This 'core' is a region where the above 
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process is not valid. The Volterra process undoubtedly introduces internal stresses in an 
elastic medium, but these stresses are non-singular on Si and S 2 insofar as d(r) obeys 
Hookean elasticity. Henceforth, the final result does not depend on the exact choice of S . 

Vortex lines: this term refers to a line defect L whose characteristic invariant is not an 
element of symmetry of the medium as above but the circulation 7 = § c u.ds of a vector u. In 
a irrotational fluid this would be the flow velocity, in a superconductor the vector potential, 

Notice that the vector u itself has the geometry of a disclination, since it rotates by a 

2n angle about the line L. One will have recognized that the singularities investigated in 
the present work are of the vortex line type, but display a geometry locally reminiscent of a 
screw dislocation in a smectic. 
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